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ABSTRACT 


Development of an effective method for evaluation of alternative energy sources in 
the automotive industry has always been a necessity for cost efficient design analysis. One 
viable alternative energy source is electricity . In the present day environment of shrinking 
fossil fuel supplies and environmental awareness, electric powered vehicles are becoming a 
low cost, non-polluting, alternative means of transportation. The analysis of reliable 
electric propulsion can be expensive without a modeling tool for evaluating design 
strategies before vehicle construction. 

This thesis explores electricity as an alternative energy source for the automobile 
of tomorrow. Under the guidelines of the World Solar Challenge, a solar powered electric 
vehicle, using a permanent-magnet brushless dc motor has been modeled and simulated in 
Simulink (Dynamic System Simulation Software). The simulations were performed with 
the goal of determining the optimum configuration to efficiently utilize the power supplied 
from the solar array, batteries, and motor. The simulated vehicle was "driven" over 
various terrain's and at various speeds. The results obtained confirm this simulation as an 
efficient design tool and present an example of an optimum vehicle speed for one 
particular vehicle configuration. 
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I. INTRODUCTION 


A. BACKGROUND 

Development of an effective method for evaluation of alternative energy sources in 
the automotive industry has always been a necessity for cost efficient design analysis. One 
viable alternative energy source is electricity . In the present day environment of shrinking 
fossil fuel supplies and environmental awareness, electric powered vehicles are becoming a 
low cost, non-polluting, alternative means of transportation. The analysis of reliable 
electric propulsion can be expensive without a modeling tool for evaluating design 
strategies before vehicle construction. 

The World Solar Challenge is a race that was designed to demonstrate current 
technology in the field of Solar Powered Electric Vehicles (SPEV). Since it's inception in 
1987, it has provided industrial and educational "teams" an opportunity to present the 
latest technological trends and breakthroughs in solar powered electric vehicle research to 
the public. The primary means of solar powered vehicle development has been 
characterized by theoretical evaluation of the vehicles velocity based on four variables: 
aerodynamic drag, rolling resistance, mass, and power output Most of the vehicles 
entered in the World Solar Challenge are evaluated with these four variables before 
construction. 

A mote effective design tool would take these variables and incorporate them into 
a dynamic vehicle model. Coupled with the motor, battery, and solar array, the design 
engineer could then base the final vehicle design on data accumulated from simulations of 
this dynamic model. A Solar Powered Electric Vehicle Simulator is therefore absolutely 
necessary as a primary design tool for cost effective electric vehicle research and 
development 

B. THESIS OVERVIEW 

The primary means for development of solar powered vehicle designs has been to 
devise a mathematical model based on the performance variables of rolling resistance, air 
drag, gravity, velocity and required system energy. This model is used to calculate the 
power required to maintain a continuous speed under steady state conditions. TTiis 
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approach is adequate for design evaluation when not considering transitions in terrain, 
speed, or cloud cover. 

The goal of this study is to develop a dynamic simulation tool for evaluation of a 
solar powered electric vehicle model based on the system variables. These system 
variables will include the vehicle's mechanical and electrical components. This approach 
will allow the designer to specify and change the vehicle's characteristics in order to 
optimize the design. In final form, this modeling tool will realistically simulate a solar 
powered electric vehicle with changing terrain, speed, and solar array power while 
providing continuous on-line data display. 

Figure 1-1 shows a block diagram of the vehicle design considered in this study. 
The primary components are solar array, battery, voltage/current bus, motor, drive train, 
and vehicle. Models for the electrical and mechanical system elements are presented 
independently and tied together to produce the system model. 

1. Background of the World Solar Challenge and Solar Powered Vehicles 

This section provides a brief history of solar powered electric vehicles and the 

events surrounding the first World Solar Challenge [Ref. 1, pp. 1-6]. Discussed are the 
basic rules governing vehicle construction and some of the vehicle designs of previous 
races. 

2. Modeling the Vehicles Electrical Components 

Discussed here are the solar array, battery, voltage/current control, and motor. 
The array and battery characteristics are developed from experimental data and placed into 
matrices which are used as look-up tables in the system model. The current bus ties 
together the array current, battery current, and the regenerative current from the motor. 
The output current form the bus is integrated to monitor the battery state of charge 
(SOC). The vehicle system model utilizes a permanent-magnet brushless dc motor which 
is presented in state space form. 

3. Modeling the Vehicles Mechanical Components 

The system variables which describe the drive train and torque loading 
characteristics of the vehicle are developed in this section. The forces on the vehicle 
consideted to be of primary importance are rolling resistance, air drag, acceleration force, 
and MU gravitational force. The coupling equations wMch describe the drive train are also 
presented in this section. 
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VEHICLE DYNAMICS 









4. Simulation Software and Model Development 

Ease of use, speed, and output capabilities are factors considered important when 
choosing the simulation software. The primary software tools were Matlab and SimuHnk 
from MathWorks [Ref. 2]. Both provide for excellent modeling of the linear and non¬ 
linear differential equations which form the basic structure of the system. Each vehicle 
component is developed and presented as a Simulink submodel. These submodels are then 
tied together to form the dynamic system of the solar powered electric vehicle. 

5. Model Response and Data Analysis 

The system model is simulated and analyzed for varied driving conditions. Plots of 
the responses are presented and discussed. The simulation results are compared to the 
theoretical expected outputs. Then, the simulation is used to present a basic driving 
strategy to optimize the chosen vehicle configuration. 

6. Conclusions and Future Work 

Finally, conclusions about the usefulness of the solar powered electric vehicle 
simulator are presented. Suggestions for expanding the model to include simulation of 
accessories, power losses, thermal characteristics, and external data accumulation 
capabilities are presented. 
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n. BACKGROUND OF THE WORLD SOLAR CHALLENGE AND SOLAR 

POWERED VEHICLES 


A. THE WORLD SOLAR CHALLENGE 

1. World Solar Challenge Race Background 

The "World Solar Challenge" race was the idea and creation of Hans Tholstrup, 
the first man to build and drive a solar powered electric vehicle 2500 miles across 
Australia [Ref. 1, p. 1-6]. He embarked upon that adventure primarily to call attention to 
the need for a reliable alternative energy source. The World Solar Challenge was an 
opportunity for private, corporate, and educational teams to prove their various vehicle 
designs and give the world an opportunity to view the technological advances being made 
in alternative energy sources for vehicles of the near future. 

The first World Solar Challenge brought a field of 24 entrants from around the 
world. It began on November 1, 1987 and was won by the General Motors team entry 
Sunracer. The average speed of the GM Sunracer was 41.8 mph, which allowed the team 
to cover the 1868 mile distance from Darwin to Adelaide in 44.9 hours. 

The design strategies used in the 1987 race set the precedent for the next World 
Solar Challenge race held in 1990, which brought 36 entrants to Australia. The course 
was again the 1868 mile Stuart Highway between Darwin and Adelaide. The race was 
won by the Biel Engineering School team entry Spirit of Biel II. The Spirit of Biel II 
averaged 40.5 mph and completed the race 46.1 hours. Table 2-1 shows the vast 
improvement in the average results achieved by all entrants who completed the race in 
1987 and 1990. 

2. Rules Governing Vehicle Design 

The general rules governing solar powered race vehicle configuration can be used 
as a basic outline for development of a modeling tool or simulation program. The vehicle 
limits considered important in this study are: 

• maximum solar collector array dimensions are 2 m wide by 4 m long by 1.6 m 
high - this limits the array power 

• maximum vehicle dimensions cannot be greater than 6 m long by 2 m wide 1.6 
m high + 1% allowance 

• minimum height of the car is 1 m 

• minimum weight of driver is 80 kg 
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These limits provide a basis for design of the submodels developed in the following 
chapters. Since solar powered electric vehicles are historically designed to be lightweight, 
the need to optimize the power delivered from the solar array is paramount. 



Average Speed, mph 

Average Finishing Time, 
h 

1987 

19 

124.4 

1990 

28 

69.8 


Table 2-1. Comparison of World Solar Challenge Race averages. 
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m. MODELING THE VEHICLES ELECTRICAL COMPONENTS 


A. THE SOLAR ARRAY COMPONENT 

The data for a solar array output to an ideal constant voltage bus was taken from 
data collection experiments performed at Rose-Hulman Institute of Technology [Ref. 3]. 
The data collection was completed by testing three types of solar array configurations; flat 
array, curved array, and flat tracking array. The following assumptions and conditions 
were present during data collection: 

• the size of the array conforms to the World Solar Challenge race rules which 
state the array must fit into a box 1.6 by 2 by 4 meters. 

• racing is done between 9 am and 6 pm. 

• the day of solar collection is clear and sunny. 

Open circuit voltages and short circuit currents were collected every 30 minutes, 
and the available power was then calculated for each array configiuation. Since the 
simulation uses the current output from an array for a specific time of day, the data may be 
programmed into a look-up table for use in the simulation routine. Figure 3-1 presents the 
normalized current output from the three array configurations versus the time of day. 
Although this data is not exact it does present an approximation that is quite close to the 
expected results as compared to the array output of the vehicles entered in the 1990 
World Solar Challenge [Ref. 1, App. 4]. Figure 3-2 shows the block diagram 
implementation of the current supplied from the solar array to the vehicle system. 

B. THE BATTERY AND BUS COMPONENTS 

The battery component of the simulation model is presented in the same form as 
the solar array. Battery discharge data supplied by the manufacturer [Ref. 1, pp. 73-75] is 
placed into a matrix and then incorporated into a look-up table for use in the simulation. 
The bus combines the solar array and battery for output to the motor, and takes 
regenerative current from the motor to the battery for charging. 

1. Modeling the Batteries and Bus 

Two types of batteries are used in the model. Lead Acid and Silver Zinc. Figure 3- 
3 shows the battery voltage versus percent charge or state of charge (SOC). The output 
voltage to the motor will have a maximum value of 100 volts, although the input to the 
voltage bus may be larger. 
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Figure 3-1. Available Current from Solar Array vs. Time of Day. 



Time in 
Hours 


Ttmeof 

DcQf 


Look-vp 

Current from 

r 

Table 

U 

Solar Array 


Bus 


Figure 3-2. Block diagram of current supplied from the solar array to the vehicle 

system. 
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The battery State-of-Charge model (SOC(t)) is developed from the Martin Model 
[Ref. 4, p. 87]. It is calculated as a function of time from the current drawn by the motor, 



SuaeofOutrge(%) 

-^LeadAcid -+-SUverSnc 


Figure 3-3. Battery Voltage vs. State of Charge. 


current supplied by the array, and current regenerated by the motor. An expression for 
this relationship is 


SOC{t) = SOC(O) - 100 


lo( ^mot ha ) 

Qo 


(3.1) 


where is defined as the battery’s maximum current capacity in amp-hours, 4^ is the 
current supplied from the solar array as defined in the previous section, and equals the 
motor current. The motor current is positive when the motor is drawing current and 
negative when the motor is generating current. This regenerative action of the motor is 
caused by slowing, stopping, or going down a hill. SOC(O) is the battery's initial state of 
charge (in percent). 
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The block diagram shown in Figure 3-4 presents the implementation of the battery and the 
voltage bus to be used in the overall vehicle model. 



Figure 3*4. Block diagram of battery and bus model to be simulated. 


C. THE PERMANENT-MAGNET BRUSHLESS DHtECT CURRENT (PMBDC) 
MACHINE 

A simplified approach to development of a permanent-magnet brushless dc 
(PMBDC) motor wiU be presented from the Krause method [Ref. 5, pp. 273-328] of 
machine analysis. This established theory of PMBDC machines wiU present the machine 
dynamic characteristics in terms of linear differential equations, which wiU then be given in 
state space form. The matrix of state equations wUl be the basis for construction of the 
time-domain block diagram for the dynamic model. 

1. The Machine Volts^e and Torque Equations 

The two-pole two-phase permanent-magnet dc machine, as shown in Figure 3-5, 
wiU be used to develop the machine equations. In the circuit diagram, represents the 
identically distributed as and bs windings on the stator and the winding resistance. The 
d axis denotes the rotor magnetic axis, and the q axis measures the angular displacement 
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of the rotor referenced from the stator as winding. The q axis is displaced 90° counter¬ 
clockwise from the d axis. The rotating magnetic field produced by the interaction of the 
stator current and the magnetic rotor develops a positive electromagnetic torque (TJ, 
which turns the rotor counter-clockwise in the direction of increasing 0^ at an electrical 
speed of co^- The load torque (T^) acts to oppose the electromagnetic torque. 

The following assumptions wiU be made in our model development: 

• No demagnetization of the permanent magnet due to large currents or 
temperature. 

• Eddy currents and harmonic losses are neglected. 

A loop analysis of the circuit shown in Figure 3-5 provides the motor voltage 
equations 


Vas=rJas+P'^<^ 

Vbs=rshs+PK 

where p is the notation used for the operator d/dt and A, is referred to as the flux linkage. 
The sinusoidal voltages and are 90° phase shifted. The flux hnkage is a function of 
the inductance of the stator windings and the magnetic flux produced by the permanent- 
magnet rotor. The flux linkage equations may be expressed as 


^as ~ ^asaJ'as ^asbs^bs ^asm, magnetic rotor 

'^bs ~ ^bsaJ'as ^bsbs^bs “^bsm, magnetic rotor 


(3.4) 

(3.5) 


and 


\_^bsm J 


= 


sin 0 ^ 
- cos 6^ 


(3.6) 


where is the amplitude of the flux linkages as seen from the stator, 0^ is defined as 


d, = Ji(0,(^)d^ + e,(0) 


(3.7) 


12 



where co^ is the electrical angular velocity and ^ is a dummy variable of integration. 

By assuming a uniform air gap and identical windings, the mutual inductance 
between the as and bs stator windings is zero and the self-inductances and are 
equal and denoted Utilizing L„ and equations (3.4) through (3.7), the voltage 
equations (3.2) and (3.3) become 

^as ^s^as P( ^ss^as ( 3 - 8 ) 

^bs = rjbs + P( LJbs + (3-9) 

or in matrix form 

V.*. = + pfL ) (3.10) 

The simulation of the permanent-magnet machine is made easier by transformation 
to a reference frame fixed in the rotor by use of the Pa±’s transformation [Ref. 5, pp. 133- 
160]. In this transformation, the electrical frequency 0)^ = (O^ and therefore 0g = 6^ Ibis 
will enable the use of a dc voltage input to control the machine speed rather than a 
sinusoidal voltage from a dc-to-ac inverter. The Park's transformation for a two-phase 
machine is 

fgsl _ 

/xj [sinQr 


( 3 . 12 ) 
and 

(3.13) 

where /represents either voltage, current, or flux linkage. The subscript s denotes stator 
variables and superscript r indicates transformation to the rotor reference frame. By 
substituting (3.13) into (3.10), the matrix voltage equation becomes 



sinQr Jfc 


cosQ^llfbs 


(3.11) 
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(3.14) 


Pre-multiplying equation (3.14) by and considering A,^ as a constant, the voltage 
equations may be written in matrix form as 


Vaj ^ "f” P^ss ^r^ss ^qs ®r^m 

+ 


v^J L ~®r4. r,+pL,,jli^ 


(3.15) 


Since 0^ = 9^ 


and equation (3.15) becomes 


(3.16) 


P^ss 


'(^r^ss PAyj Jl 


'ss ^qs 

+ 


(3.17) 


The electromagnetic torque is related to the rotor speed (O^ by 


T, =(Jr + J^)-^+ 50 ), + 


(3.18) 


and to motor current i' by 




(3.19) 


Combining equations (3.18) and (3.19), we obtain 


■“ ^m^qs = (Jr + Jml) + ^®r + ^motload 


(3.20) 
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where and represent the inertia of the rotor and the connected mechanical load 
respectively, B represents the damping coefficient due to the friction and windage losses, 
and represents the load torque delivered to the motor. 

2. The Motor as a System of State Variables 

From equations (3.17) and (3.20), expressions for the relationship between the 
voltage equations and the electromagnetic torque can be written in terms of the state 
variables i^^, ij, and co^ 





(3.21) 


dids ^ 
dt 



+ 


( 3 . 22 ) 


dm, B (pV 

dt (Jr + Jml) V2j 




1 

(Jr+Jna) 


T 

^motload 


(3.23) 


By substituting the electrical time constant Xa=L^ra, the back-emf constant 


K - torque constant k. 


th® two-phase system described 


by the linear differential equations (3.21) through (3.23), the system becomes a three- 
phase model represented by 



(3.24) 


The constants k^, and are usually provided by the manufacturer's for a specific 
motor. 

Two important observable values are the current drawn and the power delivered 
by the motor. The actual current drawn by the motor is related to the input by 
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(3.25) 


^dc^dc 



where and 1^^ are the input voltage and current to the motor. Since then the 

current drawn by the motor is given by 


4c 



(3.26) 


The actual power delivered by the motor is 


‘ motor 


^motload 


1 ? 

Pj 


CO, 


(3.27) 


The time-domain block diagram for the permanent-magnet dc machine represented 
by equation (3.24) and (3.26) is shown in Figure 3-6. The diagram represents the motor 
component of the vehicle system to be simulated. 
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Figure 3-6. Time-domain block diagram implementation of the permanent-magnet 
brushless dc motor given in equation (3.24). 
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IV. MODELING THE VEHICLES MECHANICAL CHARACTERISTICS 


A. THE FORCES ACTING ON THE VEHICLE 

The forces that have been considered important in modeling an electric vehicle are 
described as the vehicles performance variables [Ref. 7, App. B, and Ref. 1, pp. 46-47]. 
These performance variables are comprised of forces acting on the vehicle and are used to 
calculate the power required to propel a vehicle at some steady velocity. The modeling 
equation which describes this power is 

Pmotor - power required to drive the vehicle 


where, P/j = Power required at the rear drive wheel, P^ = Power required for accessories, 
Ej^ = Efficiency of the accessories, E^fEg = mechanical and electrical efficiency of the drive 
train and motor. By neglecting the accessories, we may express Pmotor ^ 


_ Pr _ v(Fw + pR Phg Pa) 


(4.2) 


where v = the vehicle velocity, F^ = wind force, Fg = rolling resistance force, FfjQ = hill 
gravitational force, and = acceleration force. These forces are describe the total drag 
on the vehicle, in Newton's and will be expanded upon. Equation (4.2) will be referred to 
as the performance equation. 

1. The Wind Force (F^) 

This performance variable describes the force of faction on the vehicle body as it 
moves through the air. It is specified at zero wind velocity and is a function of the frontal 
surface area of the vehicle, as well as any protrusions on the vehicle, air passages in the 
body and any other components of the vehicle affected by movement over, around or 
through the vehicles shell or body. The ideal vehicle shape is a long cylinder with tapered 
ends. The equation that describes the wind force is; 

j PyC,A (4.3) 
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where is the air density (1.2 kg/m^ was used in all calculations), v is the vehicle velocity 
(m/s), Q is the coefficient of aerodynamic drag, and A is the vehicles frontal surface area 
(m2). 

2. The Rolling Resistance Force iJF^ 

The rolling resistance is primarily a function of the tire friction on the roadway. 
The coefficient of friction is directly related to the type of tire used. Narrow, high 
pressure tires will have a lower coefficient of friction than tires that are wider and lower 
pressure. Bicycle tires are normally used in solar powered vehicle applications. The 
equation which describes the rolling resistance force is: 

= MyC,g (4.4) 

where My = vehicle mass (kg), = coefficient of rolling resistance, and g = gravitational 

constant (9.8 m/s2). 

3. The Gravitational Force (F^g) 

The hill gravitational force affects the total drag only when the vehicle is traversing 
a hill. Uphill equates to a positive grade and downhill a negative grade. The equation 
describing this force is: 

Frq = Myg sin 0^ (4.5) 

where 0 = angle of the road measured from the horizontal component of the driving route. 

4. The Acceleration Force (F^) 

The acceleration force is also known as the inertial force [Ref. 4, pp. 20-23]. This 
component of the total drag will have a positive sign for acceleration and negative sign for 
deceleration. The acceleration force equation is: 

F^ = Mya (4.6) 

where a = vehicle acceleration (m/s2). is a dominant force in the process of 
regenerative braking. This is the process of recovering energy in the form of current 
supplied to the battery as discussed in Chapter HI. 

5. The Load Torque to the Motor 

The forces expressed in the performance equation (4.2) relate directly to the load 
experienced by the motor. Realizing that a chain or belt drive translates these forces to 
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the motor tangent to the rotation of the drive wheel in contact with the road, the load 
torque then becomes [Ref. 8, p. 183] 



where r is the radius of the drive wheel and N is the speed or gear ratio from the drive axle 
to the motor. Expanding the sum of vehicle forces represented in equation (4.7) gives 



Figure 4-1 shows how the vehicle forces ate translated to the motor. 



Figure 4-1. Interaction of vehicle force with motor for computii^ load torque 

B. THE VEHICLE DRIVE TRAIN 

The drive train acts as the coupling mechanism between the motor and the wheels. 
The most common method used in solar powered electric vehicles are chain or belt drives. 
The drive train model will be presented by assuming a belt drive configuration. The motors 
radial velocity in revolutions per minute is given as 
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where 0)^ is the rotor electrical frequency in radians per second, and P is the number of 
poles for the specific motor. 

The RPM^f is converted to the appropriate linear velocity by multiplication with 

the gear ratio and tire radius. The selection of the gear ratio depends on the things such as 

type of terrain over which the vehicle is traveling, the desired vehicle speed, the rated 

speed of the motor, etc. For example, over a level terrain a ratio of 1:5 may be utilized. If 

on the other hand, the terrain is uphill, an appropriate gear ratio would be 1:9. A higher 

gear ratio provides more power for climbing hills, but also greatly reduces the maximum 

speed of the vehicle. Utilizing equation (4.9), the gear ratio G»=—, and the wheel radius 

N 

r, the vehicle velocity is expressed by 

Vcarikph) = RPM^, • • 2nr (4.10) 

c^rV y mot ^ 


then converting to miles per hour 


Vcar ( ) = 0. 625) 


(4.11) 


C. THE MECHANICAL COMPONENTS DESCRIBED AS A SYSTEM 

The mechanical components comprising the vehicle forces which translate into the 
load torque ^had experienced by the motor, are presented in the block diagram shown in 
Figure 4-2. The input is the angle of the incline or decline in radians. If the vehicle is 
traveling on level ground, then the input would equal zero. The vehicles forward 
velocity component of the performance equation (4.2) is derived from the drive train. The 
velocity described in equation (4.11) is converted firom kilometers per hour to meters per 
second. 
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The drive train equations (4,10) and (4,11) are presented in block diagram form in 
Figure 4-3, 



Figure 4-3. Block diagram of drive train conversion of motor rotational speed to 

vehicle forward velocity. 
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V. SIMULATION SOFTWARE AND MODEL DEVELOPMENT 


A. THE SIMULATION SOFTWARE 

Matlab [Ref. 2] was chosen as the simulation software for the vehicle system. The 
Matlab Toolbox - Simulink, is the primary building tool for designing and connecting the 
dynamic model components of the vehicle simulator. Simulink models both linear and 
non-linear system components in block diagram "windows". These windows resemble the 
block diagrams presented in many control theory textbooks, which allows the user to 
easily construct and modify a system with a click of a mouse button. The model 
components presented in the previous chapters are assembled in Simulink block diagram 
submodels, then connected to present the overall vehicle system. Figure 5-1 shows the 
Simulink model for the entire solar-powered electric vehicle simulator system. 

The user interface with the vehicle simulator is generated in MATLAB code 
(Appendix A). A user control window initiates the program and allows the user to input 
the vehicles characteristics by clicking the mouse on a "pushbutton" within the SPEV 
figure window. The required input information is displayed in the MATLAB control 
window. Once the vehicle data has been entered the simulation may be run and output 
data analyzed. The specific vehicle information required for input is discussed in Chapter 
VI. User instructions are provided in Appendix A. 

B. THE ELECTRICAL COMPONENTS IN SIMULINK 

In Chapter HI, the individual components which comprise the electrical 
characteristics of the solar powered vehicle were presented as equations and block 
diagrams. These block diagrams will now be presented in Simulink form. 

1. The Run Time and Solar Array Characteristics 

The first components of the vehicle simulator are the run time and solar array 
output. Current vectors are developed from the three types of solar array configurations 
presented in Chapter HI, which are the flat, curved, and tracking array configurations. An 
array type and expected peak current output are specified by the user, then the associated 
vector is read into the Array look-up table block. The time of day for simulation start is 
specified by the user and can be between 0900 and 1700. Figure 5-2 shows the simulink 
implementation of the time and array look-up table. The variable cloud cover block is 
used to simulate the loss of array output due to overcast conditions. The user can "double 
click" on the block, opening it to 


24 



Time in HoursN— 12:34 



Figure 5-1. Solar-Powered Electric Vehicle simulator in Simulink. 
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of important characteristics 






a slider gain. Clicking on the down or up arrow changes the gain coefficient which 
produces changes in the array current output fed to the Electrical Components block, 
thereby simulating current degradation due to cloud cover. 




Cutrentfiom 
Solar Array 

variable cloud 
cover 




hours 
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\ This block is a "Slider Gidn” 
^ch mult^lies die array current 
a coefficient between 0 and 1 


> Array curretd to ffiectricalCom|ioDaib block 


Figure 5-2. Simulink implementation of the Solar Array. 


2. The Electrical Components Block 

The Electrical Conyronents block shown in Figure 5-3 takes input from the 
vehicle throttle, which is a slider gain, current fed from the array as previously discussed, 
and current being drawn by the drive motor. The output is controlled dc voltage to the 
motor from the battery, and a multiplexed array of outputs to be monitored in the On-Line 
Display. These outputs include: 

• state-of-charge 

• total amps-hours used 

• battery voltage 

• motor current drawn 

• array current delivered 

The Battery block uses a look-up table to determine voltage output based on state- 
of-charge as presented in equation (3.1). The user specifies the type of battery to be used, 
either Lead-Acid or Silver-Zinc. The discharge curve for the chosen battery will be 
displayed in the Battery block. Also seen in Figure 5-3 are the contents of the Electrical 
Components block, which will be expanded on. 
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Figure 5-3. Simulink representation of the main electrical components of the vehicle 

simulator. 
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a. The Current Bus Block 

The Current Bus block takes input from the solar array and the current 
being drawn by the motor, then uses these to compute the battery state-of charge (SOC) 
and the amp-hours (Ahrs) used. Figure 5-4 shows the expanded view of the Current Bus 
block. As developed in chapter 3, the sum of the array current and motor current are 
integrated then, divided by 3600 seconds to provide the amp-hour usage output. Since the 
total amp-hour usage cannot be less than zero, a switch is used to "turn off' the amp-hour 
output at a value of zero. The State-of-Charge function block implements equation (3.1) 
by using the amp-hour input to compute the battery state-of-charge based on maximum 
battery capacity and initial state-of-charge. 

b. Speed and Voltage Controller Block 

The inputs to the Speed and Voltage Controller block are throttle input 
and battery voltage. The output is controlled dc voltage to the motor to regulate speed. 
Figure 5-5 shows the expanded Simulink block. The throttle input, from a slider gain is 
sent to the voltage ratio function block which outputs a coefficient from 0 to 1. The 
throttle coefficient is rate limited before being multiplied by the battery voltage. The 
mayimiim voltage to the motor is set by a saturation block at 1(X) volts. The rate liimted 
or "throttled" voltage is then sent to the motor. 

3. The Drive Motor Block 

The inputs to the Drive Motor block are controlled voltage from the Electrical 
Components block and vehicle load (Ti^^). The vehicle load will be discussed in following 
sections. The outputs are RPM to the drive train and a multiplexed array of RPM, motor 
load, efficiency and current drawn by the motor to the On-Line Display. The Drive Motor 
block is presented in Figure 5-6 and its contents are shown in Figure 5-7. 

a. The System Efficiency Block 

The System Efficiency block consists of a two dimensional look-up table 
with inputs of motor rpm and vehicle load, and switches for minimum and maximum 
efficiency. The efficiency limits are set by the minimiun and maximum efficiency for the 
specific motor being used. The System Efficiency block is presented in Figure 5-8. The 
actual vehicle load is divided by the efficiency output in the Load/Eff function block. The 
output of this fimction block relates to the T^^oad introduced in equation (3.18), and used 
in equation (3.27) to determine motor power out. This is the vehicle load torque seen by 
the motor. 
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Figure 5-4. Simulink representation of Ahrs and SOC output from current input 
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Figure 5-5. Simulink representation of Speed and Voltage control. 
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Figure 5-6. Simulink diagram of the Drive Motor block. 
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b. The PMBDC Block 

The PMBDC block represents the permanent-magnet brushless dc motor 
used in the vehicle and implements equation (3.24) in Simulink form. The Simulink 
diagram of the motor is presented in Figure 5-9. The magnitude of the dc voltage through 
input 2 directly controls the desired rpm of the motor. The actual motor speed is of 
course affected by the load torque introduced by the mechanical components of the 
vehicle and the inertial load The speed in rpm is output to the Drive Train block of 
the simulator. The motor current drawn, rpm and load torque are sent to the On-Line 
display to be monitored during simulation. 

C. THE MECHANICAL COMPONENTS IN SIMULINK 

In Chapter IV, the individual components which incorporate the mechanical 
characteristics of the solar powered vehicle were presented as equations and submodels. 
The Simulink form of these submodels will now be presented and their operation 
explained. 

1. The Vehicle Drive Train Block 

The Vehicle Drive Train block uses the rpm input to convert rotational velocity of 
the motor to the linear velocity of the vehicle. The Sim ulink system through which this is 
accomplished is presented in Figure 5-10. The rpm delivered from the Drive Motor block 
through input 1 is divided by a specified gear ratio and mechanical efficiency, then 
multiplied by the circumference of the drive wheel. This provides a linear velocity in 
meters per minute. The linear velocity is converted to meters per second, for use in the 
Vehicle Characteristics block in determination of forces, and miles per hour. The miles 
per hour variable is integrated to provide the total miles traveled. The speed, in mile s per 
hour, and total distance traveled are sent to the On-Line Display and monitored as the 
speedometer and odometer. The total miles traveled is also sent to the Vehicle 
C/iaractemftcs block for force computations. 

2. The Vehicle Characteristics (Load Torque) Block 

The Vehicle Characteristics (Load Torque) block is simply the implementation of 
the performance equation (4.2) and the associated torque equations detailed in section 4.5. 
The inputs to the Vehicle Characteristics block are the velocity in meters per second, total 
miles driven, and the voltage delivered to the motor. Each of the forces developed by the 
motion of the vehicle are summed and sent to the switch block. At this switch block the 
total force is either zero, if no voltage and no velocity are present, or the sum of all forces 


33 




Controlled 
voltage in from 
the throttle 



Current out 


Idc 


Current drawn 
by motor 




1 .. 

” ra*ta 


_^ 

l/ta 


ta.s+-l 


Electromagnetic 
Torque (Tmot) 


kt=P*ke/2 


-► 60/(P*pi) - 

■ _ ■ 

Electrical angular 
velocity 


rpml 

-*{I] 

RPMout 


— 



+ 

-^ 

?/2 

0 


H 



Jtots+B 


Vehicle load 
torque in(Tload) 


Power delivered 


* pwrdel 


Figure 5-9. Simulink representation of a Permanent-Magnet Brushless DC Motor 

(PMBDC). 
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present. Then, since the total force is actually applied at the point where the wheel 
contacts the road, the force is multiplied by the wheel radius and divided by the gear ratio, 
providing the vehicle load torque seen by the motor. The complete contents of the 
Vehicle Characteristics block are shown in Figure 5-11. Each of the force components 
within the vehicle characteristics block will be presented and expanded upon. 

a. The Rolling Resistance Block 

The contents of the Rolling Resistance block are shoAvn in Figure 5-12. 
The input to the block is the linear velocity of the vehicle and the hill force. The switches 
shown have three purposes: 1) if the vehicle velocity is zero, then the rolling resistance 
force is zero, 2) if the vehicle velocity is greater than zero, the rolling resistance force is 
applied against the direction of motion, and 3) if the hill force is greater than the rolling 
resistance force, and the vehicle is moving backwards, the rolling resistance force is 
applied in the opposite direction of motion. The rolling resistance force is a constant and 
not a function of the vehicle velocity. 

b. The Hill Gravitational Block 

The Hill Gravitational block computes the force on the vehicle as a 
function of hill angle. If the vehicle is traveling up-hill, the force is positive and increases 
the total load on the motor. If the vehicle travel is down-hill, the force is negative and 
consequently decreases the total load on the motor. If the vehicle is driving over level 
terr ain, then the hill force is zero. This component plays a significant role in the 
regenerative action of the motor. If the hill angle is large enough to overcome the other 
forces applied to the motor, then the load of the vehicle will "drive" the motor, which in 
turn wfll "generate" current for battery charging. Computation of the hiU force in 
Simulink is shown in Figure 5-13. 

c. The Acceleration Force Block 

The input to the Acceleration Force block is the linear velocity of the 
vehicle. The velocity is differentiated to obtain the vehicle acceleration used in 
determining the accelerating or decelerating force. As with the hill force, the acceleration 
force also plays a significant role in the regenerative action of the motor. If the voltage 
applied to the motor is decreased, the vehicle will begin to slow down. This action causes 
a large deceleration, or negative force, to be applied to the motor, which in turn creates a 
significant regeneration of current to the battery for charging. 
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Figure 5-12. Simulink representation of the Rolling Resistance force. 
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This force is only applied during transient operation and has no steady-state component. 
Figure 5-14 shows the Simulink representation of the acceleration force computation. 



Vehicle 

Mass 


Figure 5-14. Simulink representation of the Acceleration Force computation. 
d. The Wind Force Block 

Figure 5-15 shows the Simulink representation of the wind force. The 
input to the Wind Force block is the linear velocity of the vehicle. 



Frontal Coefficient of Density 

Area Drag 


Figure 5-15. Simulink representation of the Wind Resistance computation. 


This force is directly proportional to the square of the vehicles linear motion. The faster 
the vehicle goes, the greater the force. 
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D. THE VEHICLE ON-LINE DISPLAY IN SIMULINK 

The On-Line Display block takes inputs from various components of the 
simulation program. These inputs are all feed into a simulink multiplex block and then 
sent to a Madab s-function called spevanl.m. The function spevan2.m is included in 
Appendix B with the other Madab functions used in the vehicle simulation program. The 
On-Line Display block and its components is presented in Figure 5-16. All inputs are 
clearly labeled. 

Once the simulation routine has been initiated by the user, the On-Line Display 
window will appear at the bottom of the monitor. The outputs to the screen are clearly 
identified. Upon completion of a simulation run a "double-click" on the Plots block will 
provide a profile of important data and average values of certain characteristics. These 
average value components will be detailed in the following chapter. 
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Figure 5-16. On-Line Display block and contents in Simulink. 
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VI. TESTING AND EVALUATION OF THE SYSTEM MODEL 


A. OVERVIEW OF THE TESTING APPROACH 

The primary question to ask in testing the vehicle simulation program is, "Does the 
power out of the motor equal the power required to overcome the forces acting on the 
vehicle for a specific speed?" This chapter will be concerned with answering this question 
over three types of terrain profiles - flat, single hill, and variable elevation. 

I. Specification of the Initial Simulation Parameters 

The vehicle electrical and mechanical simulation parameters where compiled from 
actual solar powered vehicles built by teams entered in the 1990 and 1993 World Solar 
Challenge. 

The mechanical characteristics were complied from a vehicle built and raced by the 
CALSOL solar powered vehicle race team from the University of Catifomia at Berkeley. 
Table 6-1 shows the performance characteristics used for the vehicle simulator testing. 
These characteristics also serve as the default values for initiation of all simulations. 


Coefficient of Drag (Cn) 

.16 

Frontal Area (A) 

1.2 m2 

Vehicle Mass (Mv) 

295.8 kg 

Coefficient of Rolling Resistance (Cr) 

.0105 

Rear Wheel Radius (r) 

.26 m 

Gear Ratio (1:N) 

1:5 

Drive Train Efficiency 

98% 


Table 6-1. Performance characteristics used in the vehicle simulation. 

A permanent-magnet brushless dc motor built by Unique Mobility, Inc., is used in 
this vehicle simulation. The characteristics of the motor used in the simulation are 
presented in Table 6-2. Figure 6-1 shows the efficiency map provided for the motor. This 
efficiency map is loaded into the program upon initi ali zation of the motor parameters for 
use during simulation. As with the vehicle characteristics, the motor parameters used for 
simulation testing also serve as the default values for simulator initialization. 
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Unique Mobility DR127s/CR10-100 

Voltage In 

50 volts 

Rated Power 


Max No-Load Speed 


Number of Poles (P) 

24 

Damping Coefficient (B) 

0 

Rotor Inertia (J) 


Electrical Time-Constant (T^) 

6.944e-4 

Winding Resistance (ra) 

3.6e-2 ohms 

Back-EMF Constant (k^) 

20 V/Krpm 

Torque Constant (k^) 

.191Nm/A 


Table 6-2. PMBDC Motor parameters used for simulation testing. 


Efficiency plot for Unique DR127/10 



Figure 6-1. Motor efficiency plot 
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The battery chosen for simulator testing is a Lead-Acid type, with a capacity of 35 
Ahrs and initial state-of-charge of 100%. A tracking solar array will be used, with a peak 
current output of 11 amps. This battery and array configuration are the default 
components used in simulation initialization. 

B. MODEL SIMULATION 

The three types of terrain considered during testing are flat, single hill, and 
variable terrain which contains numerous positive and negative grades. The following 
outputs will be presented: 

• driving terrain profile 

• vehicle speed 

• sum of the forces on the vehicle 

• current drawn by the motor 

• torque seen by the motor 

• power delivered by the motor 

• power required at the rear drive wheel 
1. Simulation Over Flat Terrain 

The first driving profile for simulation is over flat terrain as shown in Figure 6-2. 
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Figure 6-2. Plot of flat terrain driving profile. 

The simulator system was initialized with the parameters specified in the previous section. 
Cloud cover was neglected, so peak current from the solar array was 11 amps. 
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The first outputs examined are the vehicle speed and total force on the vehicle. 
The theoretical speed is calculated from a quadratic equation based on the relationship 
between the energy delivered by the motor and the energy required to overcome the 
vehicle forces. This is expressed as 


A + B(rprnl^t) = C - D(rpm^,) 

_ MygCff _ p^CjyA2n^r^ 

NE^ ’ 60 ^N^E^ 


( 6 . 1 ) 



With data provided in Tables 6-1 and 6-2, we get 


+ 


0.1055 


0. 000000607 




^ (6.2- 263.9) _ ^ 
0. 000000697 


Then using the quadratic equation 


( 6 . 2 ) 




_ -b ± - 4ac 

2a 

= -153770 or 2404. 5 


(6.3) 


The motor rpm in (6.3) equates to a vehicle speed of = 29.4 mph or 13.1 m/s. lliis 
value is within 95% of the simulated value presented in Figure 6-3. The slight variation is 
most likely caused from not considering the motor efficiency in the calculation. 

Using a speed of 13.1 m/s, the theoretical total steady state force acting on the 
vehicle is derived from equation (4.2) as 


^vehicle = K + P'r 

= MygCji -t- — p^CjjAv^ 

= (295.8)(9. 8)(0. 0105) + -(1. 2)(0.16)(1. 2)(13.1^) 

2 

= 50. 2 ( Newtons ) 

(6.4) 
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Figure 6-3. Plots of the steady state speed and total force acting on the vehicle for 

flat terrain driving 
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The acceleration force (F^) and hill gravitational force (F^^) are zero for flat steady state 
driving. The force given by (6.4) is within 96% of the simulated steady-state value of the 
total force displayed in Figure 6-3. 

Also observed during the simulation are the current drawn by the motor and the 
torque delivered to the motor. Using equation (4.8), the theoretical force provided in 
(6.4), and data given in Table 6-1, the theoretical load torque delivered to the motor is 



{0. 26)(50,2) 
(5)(0.98) 


= 2.6(Newton - meters) 


(6.5) 


Utilizing the motor rpm from (6.3) and the load torque from (6.5), the motor efficiency 
may be determined from the efficiency plot in Figiue 6-1. The efficiency is approximately 
75%. Dividing (6.5) by the efficiency, the load torque delivered to the motor is 


0.75 


= 3. 47(Newton - meters) 


Then utilizing (6.6) and equation (3.24) the theoretical current drawn by the motor is 


As seen in Figure 6-4, the theoretical load torque and motor current drawn correspond 
closely to the steady-state output from the simulated model. Both values are within 97% 
of expected results. 

The most important comparison to make is the power required to drive the vehicle 
versus the power delivered by the motor. The theoretical power required is computed 
from the performance equation (4.2). Utilizing the vehicle velocity from (6.3) and the 
total force on the vehicle given in (6.4), the power required at the rear wheels is 

P^=v(F^ + F^) = (13.1)(50.2) = 657.62( Watts) (6.8) 


Dividing (6.8) by the motor efficiency and the drive train efficiency, the power 
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Figure 6-4. Plot of the motor current and the load torque delivered to the motor for 

flat terrain driving. 
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required from the motor is 


R _ _ 




657. 62 
(0.75)(0.98) 


= 894. 72 W 


The output plots of the power delivered by the motor and the power required at 
the rear wheels with consideration for system efficiencies are presented in Figure 6-5. The 
simulated steady-state values are within 94% of the computed power in (6.9). 

2. Simulation Driving Over a Single Hill Profile 

The next simulation was up a single hiU. Again, the array current was a constant 
11 amps. The terrain profile for this simulation is shown in Figure 6-6. Using (6.1), with 


^ ^ MygCg ^ MygsinQ^ii 


( 6 . 10 ) 


the speed for the terrain profile with 0^; = 0. 03789 radians is 


2 0.1055 (28.62- 263.9) _ 

+ Q 000000697 0. 000000697 


( 6 . 11 ) 



( 6 . 12 ) 


The motor rpm in (6.12) equates to a vehicle speed of = 2<5,9 mph or 12.0 m/s. This 
value is within 94% of the simulated value presented in Figure 6-7. 

Using a speed of 12.0 m/s, the theoretical total steady state force acting on the 
vehicle is 

Vehicle ~ ^R ^HG ~ ^vS^'R + ^v8 ^hUl 

= (295. 8)(9.8)[(0. 0105) + sin(0. 03789)] + | f 2)( 0.16)(1. 2)(12. 0^ ) 

= 156. 8(Newtons) 

(6.13) 
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Figure 6-5. Plots of the power required at the rear wheel and power delivered hy 

the motor for flat terrain driving. 
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Figure 6-6. Plot of single hill terrain profile. 
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Figure 6-7. Plots of the vehicle speed and the total force acting on the vehicle for a 

single hill terrain profile. 
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The simulated force in Figure 6-7 is within 97% of the computed force in (6.13). 

Using equation (4.8), the force from (6.13), and data from Table 6-1, the 
theoretical torque delivered to the motor for this terrain profile is 


_ r(Fw + pR + ^HG) 

^ (0. 26)(156.8) 
(5)(0.98) 

= 8. 32(Newtons - meters) 


(6.14) 


I T rill zin g (6.12), (6.14), and the efficiency plot in Figure 6-1, the motor efficiency is 
approximately 85 %. By dividing (6.14) by the motor efficiency, the theoretical load 
torque delivered to the motor is 

Tnu>tioad = 7^ = 9.78(Newton - meters) (6.15) 

U, oj 


Then utilizing (6.15) and equation (3.24) the theoretical current drawn by the motor is 



51. 2( Amps ) 


(6.16) 


As seen in Figure 6-8, the theoretical load torque and current are within 92% of the values 
obtained from the simulation. 

For the power comparison, the theoretical power required for this driving profile is 
computed from the performance equation (4.2). Utilizing the total force on the vehicle 
given in (6.13) and the vehicle velocity, the power required at the rear wheels is 

Pg=v(F„ + Fg + Fg^) = (12.0)(156.8) = 1881.6(Watts) (6.17) 

Dividing (6.17) by the motor efficiency and drive train efficiency, the power required from 
the motor is 


_ Eg 


1881.6 

(0.86)(0.98) 


= 2258.5( Watts) 


(6.18) 
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current drawn by motor (A) 



Load torque in (Nm) 



miles driven 


Figure 6-8. Plots of the current drawn by the motor and the load torque delivered 
to the motor for a single hill driving profile. 
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The output plots of the power delivered by the motor and the power required at the rear 
wheels with consideration for system efficiencies are presented in Figure 6-9. The 
simulated values are within 92% of the expected value of power computed in (6.18). If 
the vehicle speed had been computed with consideration of motor efficiency, there would 
have been less variation between theoretical and simulated results. 

3. Simulation Driving Over Variable Elevation Terrain 

The next simulation is a variable or "hiUy" terrain. Solar array current was a 
constant 11 amps. The terrain profile for this simulation is shown in Figure 6-10. 

Since the simulation has been shown to correspond closely with expected values 
for the flat and single hill driving profiles, the plots for this variable terrain will be 
presented without theoretical comparison. As seen in Figures 6-11 through 6-13, the 
outputs are consistent with the results comparisons made in the two previous driving 
profiles. Observation of Figure 6-13 shows that the power delivered equals the power 
required. This would indicate that the simulation model does accurately portray an 
electric vehicle driving over a variable terrain. 

Also observed in Figure 6-13 is the regenerative action of the motor when the 
forces acting against the vehicle motion are overcome by forces which act to propel the 
vehicle forward. This creates a current regeneration action from the motor which is used 
by the system to charge the battery, increasing the battery state-of-charge. 

4. The Optimum Balance Between Motor Current and Speed 

Since the objective of this simulation model is to optimize the design of the 
vehicle, the next test should be to find the most efficient relationship between current 
drawn by the motor and speed of the vehicle. This test was driven over flat terrain and the 
voltage to the motor was stepped in increments of 10 volts using the throttle control. At 
each voltage, the system was allowed to settle to steady state, then the current being 
drawn by the motor and the vehicle speed were recorded. Table 6-3 shows the data 
accumulated from this experiment 

Since the current drawn by the motor is a function of the motor efficiency, there 
should be a range of speed in which the motor is most efficiently using the energy 
supplied. By plotting the current as a function of throttle voltage and speed, the most 
efficient operating points may be easily identified. By examining the output relationship 
shown in Figure 6-14, the optimum driving speed for this particular vehicle configuration 
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power power 


Power required at the rear wheel 
considering motor efficiency (W) 



miles driven 


Power delivered to the rear wheels (W) 



miles driven 


Figure 6-9. Plots of the power required at the rear wheel and the power delivered 
by the motor for single hill driving profile 
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Figure 6-11. Plots of the vehicle speed and the total force acting on the vehicle for a 

variable terrain profile. 







current drawn by motor (A) 



Load torque (Nm) 



Figure 6-12. Plots of the current drawn by the motor and the load torque delivered 
to the motor for a variable terrain profile. 
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power power 


Power required at the rear wheel 
considering motor efficiency (W) 



Power delivered to the rear wheel (W) 



Figure 6-13. Plots of (he power required at (he rear wheel and the power delivered 
by the motor for a variable terrain profile. 
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would be between 25 and 30 mph. This test could also be used to determine the optimum 
driving speeds for variable terrain's and different gear ratios. 


Throttled Voltage (V) 

Motor Current(A) 

Vehicle Speed (mph) 

14 

26.5 

7.5 

24 

21.4 

13.5 

34.2 

18.9 

19 

44.2 

17.4 

24.9 

54.3 

17.35 

30.2 

64.4 

18.5 

35.5 

74.5 

20.3 

40.1 

84.5 

22.6 

44.9 

94.7 

24.9 

49.1 

100 

26.1 

51.2 


Table 6-3. Motor Current and Vehicle Speed for Incremental Voltages. 

Efficient vehicle configurations can be achieved by changing the individual vehicle 
parameters necessary to optimize the current - speed relationship. 
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vn. CONCLUSIONS 


There is a need for efficient analysis of alternative energy propulsion systems for 
the automobiles of tomorrow. By developing systems that quickly analyze a particular 
vehicle configuration, the system designer can save the time and money required to 
actually build full scale models. This thesis has provided a method for the efQcient 
evaluation of electricity as an alternative energy source in automobile propulsion. 

A. CURRENT APPLICATION 

The original goal of this thesis was to model a solar powered electric vehicle for 
use in design analysis of vehicles configured under the rules of the World Solar Challenge. 
That task has been successfully completed. 

Presentation of this thesis was made to the Solar Phantom Team at Rose-Hulman 
Institute of Technology in Terra Huete, Indiana on March 6, 1995. The system evaluation 
and output presented in the previous chapters were verified by observations made by the 
Rose-Hulman team during previous solar powered races in the Solar Phantom I and 11. 
The team is presently constructing the Solar Phantom III, which has a radically different 
design than the previous two cars. The simulation model developed in this thesis wiU be 
used to develop a daily driving strategy for the Solar Phantom Team when they participate 
in Sunrayce 95 on June 20, 1995. This strategy should assist the team in efficiently using 
the energy available from the solar array, battery and motor at variable speeds and over 
variable terrain. 

B. FUTURE WORK 

This system is by no means complete. Future work could include improvement in 
the areas of: 

• modeling temperature effects on the motor and solar cells 

• modeling the motor to account for harmonic losses 

• modeling a friction braking mechanism 

• adding regenerative current control 

• incorporating the effects of auxiliary components in power consumption 

• incorporate a random number generator for the cloud cover coefficient to 
simulate variability of sunlight to the solar array 

Other work which could make this system more ^plicable to industry would be to 
develop the models for major components such as the battery, motor, vehicle 
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characteristics, and solar array as individual "plants" and evaluate them with a digital 
Signal Processing and Control Engineering (dSPACE) system. The plant models could be 
built in actual hardware, then connected to the simulation program through the d-Space 
interface. Through this connection, the simulation program may be used to operate the 
hardware plant. The systems designer may observe how hardware being considered for 
use in the vehicle design, reacts under actual operating conditions. The d-Space model 
would be much more cost efficient for actual component evaluation than construction of 
full scale models. 
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APPENDIX A. SIMULATOR USER INSTRUCTIONS 


1. The Solar/Electric Vehicle Simulator will only run from a Madab control window. The 
specific version of Madab is at least 4.1 with Simulink version 1.3c. Ensure that these are 
loaded on your system before attempting to run the simulator. 

2. Once in the Madab control window, type - solar. The user control window and the 
default terrain profile figures should appear in the upper right comer of the monitor. 
Figure A-1 shows the user control window for the simulator. 


Solar/Electiic Vehicle Simulator 


FQe Edit Window Help 





Ijjgg’lTMlI 

1 w \ •;,] 1 





1 










Figure A-1. SPEV simulation window 


Figme A-2 shows the default driving terrain profile loaded upon simulator system 
initialization. Each button in the user control window invokes a specific function. Theses 
functions are explained below 

a) Load finitial Data: Loads the initial default values for all electrical and 
mechanical parameters. Once loaded, aU values are displayed in the Madab control 
window. 

b) Run Time: Modifies the total simulation time, time of day to start driving, and 
integration time step. 

c) Motor Parameters: Allows the user to modify or change the motor 
characteristics, or implement a new motor. Input the requested data at the prompts in the 
Madab control window. 
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d) Vehicle Parameters: Allows any of the vehicles mechanical characteristics to 
be modified. Input requested data at the prompts in the Madab control window. 


- Tenain Profile 

File Edit ^^iidow He)^ 



Figure A-2. Default driving terrain profile 


e) Array Type: Change between tracking, flat, or curved solar array 
configurations and specify the peak array current 

f) Battery Type: Change between Lead-Acid or Silver-Zinc battery packs. 

g) Terrain Profile: Allows the user to choose one of three terrain profiles. At 
the prompt in the Madab control window choose flat single hiU, or load a file of variable 
terrain data. If the variable terrain is chosen, the data must be entered from a 
"filename.dat" text file with two colunms of information. Colunm one has distance in 
increments of .1 miles and column two contains the corresponding elevation in feet above 
sea level Each time a terrain profile is changed, a modified terrain profile window is 
displayed. 

h) Show All Data: Provides an updated display of all important parameters 
loaded into the simulation program. This function useful in verifying changes in system 
parameters. 

i) Motor Efficiency Plot: Provides a 3-dimensional plot of the motor efficiency 
as provided by the manufacturer. The user must supply the efficiency matrix for the 
specified motor. The columns correspond to increasing load torque, while the rows 
correspond to increasing motor rpm. 

j) Welcome: Shows a brief welcome to the simulator. 
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k) Help: Provides a brief help on running a simulation from the simulation 
window. 

l) Close: Closes the user window. 

m) Open Simulation Window; Opens the primary vehicle simulation window. 
Once all previous parameter input requirements have been meet, simply click on the 
"Simulation - Start" button. 

3. Once the simulation has begun, the data output window shown in figure A-3 will be 
displayed at the bottom left of the monitor. 


— 

Output 

Fi] 

le Edit 

V^ndow Help 




Motor Current (amps): 

• 

Aim Used: 

Array Current (amps): WM: 

Run time (secs): Hill Grade (%); 

Odometer (miles): 

Speedometer (mph): 

System Efficiency (%): 


Mph 

^ 30 

SOC(%); 

Battery Vohage (vohs): 
Vehicle Load (N-m): 

\ 

45 


Figure A-3. Simulation On-Line Display window 


The output indicators are self-explanatory and wiU provide continuously updated data 
display throughout the simulation. 

4. The simulation may be stopped at any time simply by clicking on the "Simulation - 
Stop" button in the main Simulink simulation window. 

5. By double-clicking on the "Plots" box in the Simulink simulation window, plots of the 
important vehicle characteristics, as well as average values of output data wiU be displayed 
in four separate figure windows. Each of these windows may be closed by clicking on the 
"File - Close" indicator. 
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APPENDIX B. PROGRAM FUNCTIONS 


A. FUNCTION SOLAR.M FOR MAIN PROGRAM 

% Solar Powered Electric Vehicle Simulator 
% Masters Thesis EE Naval-Post Graduate School, Monterey 
% This system will simulate a solar/electric powered vehicle 
% for any type of vehicle mechanical or electric characteristics. 
% 

% This system may be used to optimize vehicle configurations. 
% each function is explained separately... 

% 


% This function initiates the system user control window and 
% the default terrain profile. 

% 


% This system was developed by Steven Roerig, LT, USN 
% Date: 5 Mar 1995....all rights reserved... 


% 


% Create the Efficiency matrix (call efffunc) 

% initialize the all data matrices and set default terrain 
% to flat. 


res=10; % 10 times size of efficiency matrix 

nmot = 0; % initialize matrix 

[effmat rpmin Tlin nmot isnew]=efffunc(res,nmot); 


datamat = zeros(3,23); 

% parameter initialization 

strmat = [' 


hildef=0; 



[distvec hiUang initelev hUlvec hildef]=hillfunc(hildef); 




% Initialize the user control window and set the position 
% 

fignumber=figure(... 

’Position',[0 0 570 200],... 

’Color',[001],... 

’NumberTide','off,... 

’Name','Welcome to the Solar/Electric Vehicle Simulator’,... 
'Resize','off,... 

’Pointer','arrow'); 


axes('Units','normalized',... 
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•Position',[.15 .25 .75 .45]); 

plot([0 2 2.2.2.4 2.4],[0 0.1.1.2 .2],'m',... 

[3 3.4 5.4 5.2 3.2],[-.01.2 .2 .1 .l],'m’,... 

[8 6 6.2 7.4],[0 0.1 .l],’m’,[6.2 6.4 8.4],[.1.2.2],'m', 
[9.4 9 11],[.21-.01.22],'m',... 

'LineWidth',12) 
axis([-.5 11.5 -.1.3]) 
axis off 


^ ******************** He * ********** * 

% Create the OPEN SIMULATION WINDOW Button 


spevbut = uicontroK... 
’Style',’pushbutton',... 

'Units',’pixels',... 

'Position',[320 10 16030],... 

'String','Open Simulation Window',... 
•Visible','off,... 

Enable','off,... 

'Callback','newsim'); 


9{c % 9fe sfc 9{c »|c 9fc 9{e ^ 3{e * He ******************** * 

% Create WELCOME Button 


titll = Welcome'; 

welcscm =[' * Welcome to the Solar Vehicle Simulator. 
' * Please make your selections from the 
' * Control Buttons in the window and follow ’ 

' * the directions at each prompt in the ' 

' * Matlab Command window. ' 

» * * 

» * * 

» * * 

' * -- Enjoy the ride! - ']; 


welcbut = uicontrol(... 

'Style','pushbutton',... 
'Units','pixels',... 

'Position',[5 10 100 30],... 
Visible’,'off,... 

'String',Welcome',... 

Enable','off,... 

'Callback','helpsp(titl 1 ,welcscm)'); 


^********** ************************* 
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% Create the HELP Button 
% 

titl2 = 'Solar/Electric Vehicle Simulator Help'; 
helpscm = [' ’ 

' The SPEV Simulator will assist you in ' 

' optimzing the configuration of your ' 

' vehicle by providing user controlled ' 

' parameter variation and continuous ' 

' display. Your vehicle characteristics ' 

' may be easily changed without the initial' 

’ design costs. ' 

I t 

' Use the "User Control Buttons" to input ' 

' the vehicle characteristics, then change' 

' individual characteristics by following ' 

' the directions in the Matlab Control ' 

' Window. 

I I 

' Click on "Open Simulation Window" to open' 

' the Simulink model for the vehicle. ' 

I I 

' Use the "Throttle" to adjust the vehicle' 

' speed and "Cloud Cover" to adjust the ' 

' current from the solar array (simulate ' 

' cloud cover)... see the results. ']; 

helpbut = uicontrolf... 

'Style','pushbutton',... 

'Units'j'pixels',... 

’Position',[110 10 100 30],... 

Visible','off,... 

'String','Help',... 

'Enable','off,... 

'Callback','helpsp(titl2,helpscm)'); 

9|e 9|e s)e :(c 9|c % 4:3(e :fe 9fe 9|c He 3fe sje He He 3^ 3ie sK ^ ^ sfc 9ic 

% Create the CLOSE Button... closes the user control window 

closebut = uicontrolf... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[215 10 100 30],... 

Visible','off,... 

'String','Close',... 

Enable','off,... 
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'CaUbackVclose(gcf)'); 


^H«****H«**si« *************** ********** 


% Create all User Control Buttons 


% LOAD INITIAL DATA... loads all initial and default data for 
% first run simulations. 

% The default is mechanical data from UC Berkeley’s solar car 
% with a Unique Mobility 11.3 Hp motor, tracking array, and 
% lead-acid batteries. 

% 

firstbut = uicontrol(... 

'Style','pushbutton',... 

'Units','pixels’,... 

’Position',[5 170 140 30],... 

Visible','off,... 

'String’,'Load Initial Data',... 

Enable','off,... 

'Callback','[datamat strmat]=spevfunc("first",datamat,strmat);'); 


% RUN TIME... controls the total run time desired, the time 
% of day to start simulation, and the integration time step (dt) 

% which controls the speed of simulation. 

% 

runbut = uicontrol(... 

'Style','pushbutton',... 

'Units','pixels',... 

’Position’,[150 170 100 30],... 

'Visible’,’off,... 

'String',’Run Time',... 

Enable','off,... 

'CaUback','[datamatstrmat]=spevfunc("runtime",datamat,strmat);'); 


% MOTOR PARAMETERS... allows the user to change the motor type or 
% individual parameters ... default is a 11.3Hp/8.4KW motor. 

% 

motbut = uicontrol(... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[255 170 145 30],... 

'Visible','off,... 

'String','Motor Parameters',... 
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'Enable','off,... 

'Callback','[datamat strmat] =spevfunc("motor",datamat,straiat);'); 


% VEfflCLE PARAMETERS ... allows the user to change the vehicles 
% mechanical parameters. The default data is obtained from the 
% UC Berkeley (CALSOL) team vehcile. 

% 

vehbut = uicontrol(... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[405 170 155 30],... 

Visible','off,... 

'String','Vehicle parameters',... 

Enable','off,... 

'Callback','[datamatstrmat]=spevfunc("vehicle",datamat,strmat);'); 


% ARRAY TYPE... allows the user to choose between tracking, 
% flat, or curved array configurations and specify the peak current 
% expected. Default is tracking @11 amps peak 

% 

araybut = uicontrol(... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[5 135 140 30],... 

'Visible','off,... 

'String','Array Type',... 

Enable','off,... 

'Callback', '[datamat strmat]=spevfunc("saray",datamat,strmat);'); 


% BATTERY TYPE... choose between lead-acid or silver-zinc batteries. 
% battery voltage is set to maintain a 100 volt bus voltage. 

% Default is lead-acid battery. 

% 

battbut = uicontroK... 

'Style','pushbutton',... 

'Units','pixels',... 

Eosition',[150 135 100 30],... 

Visible','off,... 

'String','Batt Type',... 

Enable','off,... 

'Callback','[datamatstrmat]=spevfunc("battery",datamat,strmat);'); 
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% CHANGE TERRAIN PROFILE... allows user the change the driving terrain. 
% Choose between flat, single hill, or terrain profile from a data file. Data 
% must be loaded from a 2 column "filename.dat" text file. First column 
% is distance from start in 1/10 mile increments and the corresponding 
% elevation is feet above sea level. Default is flat terrain. 

% 

hillbut = uicontroK... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[255 135 145 30],... 

Visible','off,... 

'String',' Change Terrain Profile',... 

'Enable','off,... 

'Callback','[distvec hdlang initelev hiUvec hildefj=hillfunc(hildef);'); 


% SHOW ALL DATA... shows all significant data into simulator. 
% Useful after making changes. 

% 

showbut = uicontroK... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[405 135 155 30],... 

'Visible','off,... 

'Strmg’,'Show All Data',... 

'Enable','off,... 

'Callback','[datamat strmat]=spevfunc("showall",datamat,strmat);'); 


% MOTOR EFFICIENCY PLOT... 3-dimensional plot of motor efficiency. 
% Each motor must have an efficiency matrix. 

% 

effbut = uicontroK... 

'Style',’pushbutton',... 

'Units','pixels',... 

'Position',[405 100 155 30],... 

'Visible','off,... 

'String',’Motor Efficiency Plot',... 

Enable','off,... 

'CaIlback','e^ltf(rpmin,Tlin,effmat,isnew);'); 


% CHANGE MOTOR EFFICIENCY... change motor efficiency. 
% Each motor must have an efficiency matrix. 
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% 

effmotbut = uicontrol(... 

'Style','pushbutton',... 

'Units','pixels',... 

'Position',[235 100 165 30],... 

'Visible','off,... 

'String','Change Motor Efficiency',... 

'Enable','off,... 

'Callback','[effmat rpmin Hin nmot isnew]=efffunc(res,nmot);'); 

% Turns on all buttons .... 

% 

set([welcbut spevbut helpbut closebut firstbut runbut motbut... 
vehbut araybut battbut hillbut showbut effbut effmotbut],... 
'Enable',’on',... 

'Visible','on’); 


B. FUNCTION SPEVFUNC.M TO CHANGE SYSTEM PARAMATERS 
% Solar Powered Electric Vehicle Simulator. 

% Developed by Steven J. Roerig, LT, USN, at the Naval-Post Graduate 
% School, Monterey, Ca. Date: 7/26/94 
% 

% This function will change the system paramters when called. 

% 


function [spevdata,strdata] = spevfunc(action,data,strmat); 


% Load the initial data and vehicle parameters... 

% 

if strcmp(action,'first'); 

dispCSince this is your initial run, the default data has been entered') 
dispCfor a 11.3 Hp, tracking array, Pb-Acid vehicle'); 
data(l,(l;10)) = [.5 3600*9 9 24 0.014.036 25e-6/.036 20.184]; 
data( 1,(11:23)) = [.16.0105 1.2 295.8 .98 35 100 0 0.26 5 0 0]; 
data(2,:) = 1.1*[10 10 10 10 10 10 10 10 10 10... 

10 10 10 10 10 10 10 9.9 9.8 0 0 0 0]; 
strmat(l,:) = 'Tracking '; 

data(3,:) = 60*([.6.9 1.25 1.4 1.55 1.65 1.7 1.75 1.8 1.84 1.88 ... 

1.91 1.94 1.96 1.98 2.0 2.02 2.04 2.06 2.07 2.08 0 0]); 
strmat(2,:) = 'Pb-Acid ’; 

end 


% Modify the run time parameters. 
% 
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if strcmp(action,'runtime') 

dispCNote: All run times will begin at 9 AM. Maximum run time') 
dispC'is 8 Hours.') 

tmin = inputC Enter length of simulation run (minutes):'); 
data(l,2) = 60*tmin; %final time in seconds 

timein = inputC Simulation start time (btwn 0900 and 1800):'); 
data(l,3) = timein/100; 

data(l,l) = inputC Time step (how slow/fast to run: btwn .0001 and 1 sec):'); 
end 


% Modify the motor characteristics 
% 

if strcmp(action,'motor') 
dispCMotor Characteristics -') 
data(l,4) = input(' Number of poles (P):'); 
data(l,5) = inputC Damping coefficient (B):'); 
data(l,6) = input(' Inertia (J):'); 
data(l,7) = inputC Winding resistance (ra):'); 
data(l,8) = inputC Electrical time-constant (ta):'); 
data(l,9) = inputC Back-EMF constant (ke):'); 
data(l,10) = inputC Torque constant (kt):'); 
end 


% Modify the vehicle characteristics 
% 

if strcmp(action,’vehicle') 
dispCVehicle Characteristics -') 
data(l,l 1) = inputC Drag coefficient (Cd):'); 
data(l,12) = inputC Rolling resistance coefficient (Cr): *); 
data( 1,13) = inputC Frontal area (A in m'^2):'); 
data(l, 14) = inputC Vehicle mass (Mv in kg):'); 
data(l,20) = inputC Drive wheel radius (r in meters):'); 
data(l,21) = inputC Drive ratio (belt, chain, or gear - 1:N) (N=):'); 
effm = inputC Drive train efficiency (%):'); 
data(l,15) = effm/100; 
end 


% Change solar array type and peak output current 
% 

arr.dat = [10 10 10 10 10 10 10 10 10 10... 

10 10 10 10 10 10 10 9.9 9.8 0 0 0 0]; 
if strcmp(action,'saray') 
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ar_choiee = input('Solar Array type: Flat(l),Curved(2),Tracking(3): Number -'); 
peak = input('What is your expected peak current:'); 
if ar_choice = 1 

data(2,:) = (peak/10)*[.72 .78 .83 .88 .92 .97 .98 .98 .99 .985 .98 ... 

.95 .92 .90.85 .82.77 .67 .48 0 0 0 0].*arr_dat; 
strmat(l,:) = 'Flat 
elseif ar_choice == 2 

data(2,:) = (peak/10)*[.59 .63 .68 .75 .79 .85 .87 .87 .88 .87 ... 

.86 .85 .82 .78 .72 .69 .64.55 .48 0 0 0 0].*arr_dat; 
strmat(l,:) = 'Curved 
elseif ar_choice == 3 
data(2,:) = (peak/10)*arr_dat; 
strmat(l,:) = 'Tracking '; 

end 


end 


% Change battery type. 

% 

if strcmp(action,'battery') 

bat_choice = input('Battery type: Pb-Acid(l),Ag-Zn(2): Number -'); 
if bat_choice == 1 

data(3,:) = 60*([.6.9 1.25 1.4 1.55 1.65 1.7 1.75 1.8 1.84 1.88 ... 
1.91 1.94 1.96 1.98 2.0 2.02 2.04 2.06 2.07 2.08 0 0]); 
strmat(2,:) = 'Pb-Acid '; 

elseif bat_choice == 2 

data(3,:) = 68*([.8 1.35 1.45 1.47 1.48 1.49 1.495 1.497 ... 

1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.55 1.7 1.8 1.85 1.9 0 0]); 
strmat(2,:) = 'Ag-Zn 
end 

data(l,16) = inputCWhat is the A-h rating of your battery:'); 
data(l,17) = inputCWhat is your initial state of charge (percent):'); 
end 


% Show all data upon initialization and whenever called. 

% 

if strcmp(action,'showall') I strcmp(action,'first') 
dispC) 

dispCThe below information is currently is loaded into the simulator:') 
disp([' Run time: ',num2str(data(l,2)/60),' minutes']) 
disp(['Start time: ',num2str(data(l,3))]) 
disp([' Time step: ',num2str(data(l,l))]) 


dispC) 
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dispCMotor Characteristics') 
disp([' Number of poles: ',num2str(data(l,4))]); 
disp([' Damping coefficient: ',num2str(data(l,5))]); 
disp([' Inertia: ',num2str(data(l,6))]); 
disp([' Armature resistance: ',num2str(data(l,7))]); 
disp([' Electrical time-constant: ',num2str(data(l,8))]); 
disp([' Back-EMF constant: ',num2str(data(l,9))]); 
disp([' Torque constant: ',num2str(data(l,10))]); 

dispO; 

dispCVehicle Charcteristics'); 

disp([' Drag coefficient: ',num2str(data(l,l 1))]); 

disp([' Rolling resistance coefficient: ’,num2str(data(l,12))]); 

disp([' Frontal area (meters sqrd): ',num2str(data(l,13))]); 

disp([' Vehicle mass (kilo-gram): ',num2str(data(l,14))]); 

disp([' Drive wheel radius (meters): ',num2str(data(l,20))]) 

disp([’ Drive ratio (N): l:',num2str(data(l,21))]); 

disp([' Drive train efficiency: ',num2str(data(l,15)*100),'%']); 

dispC’); 

disp(['Solar Array Type: ’,strmat(l,l:21)]); 

dispO; 

disp(['Battery Type: ',strmat(2,l:21)]); 


dispC) 


dispCIf you need to make changes in your data, just click on the') 
dispCitem in the Menu Bar.') 
end 


% Output to the system. 

% 

spevdata = [data(l,l:23);data(2,l:23);data(3,l:23)]; 
strdata = [stnnat(l,l:23);strmat(2,l:23)]; 

% end of function 


C. FUNCnON HILLFUNC.M TO CHANGE TERRAIN PROFILE 

% This function performs all hill and terrain profile modifications 
% when called from the user control window. 

% 


function [distvecjiillang,initelev,hillvec,hildef]=hillfunc(hildef); 
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ifhildef~=0 

whichone = inputCFor a single hill enter -1, hill profile from file enter - 2, or flat driving 
enter - 3:'); 

dispC) 

else 

whichone = 3; 
end 


% For a single hill terrain profile. 

% 

if whichone == 1 

stl = input(’How many miles is the hill from simulation start?'); 

st2 = inputCHow long is the hill (in miles)?'); 

ell = inputCWhat is your initial evelation (feet)?'); 

el2 = inputCWhat is will be your final elevation?'); 

distvec = 0:.l;stl+st2; 

if(el2-ell)>0 

hillvec = [ell*ones(size(l:length(0:.l:(stl-.l)))) ell:((el2-ell)/(10*st2)):el2]; 
elseif (el2-ell) <0 

hiUvec = [ell*ones(size(l:length(0;.l:(stl-.l)))) ell:-((el2-ell)/(10*st2)):el2]; 
end 

type = 'Single Hill'; 

% or for a flat driving profile 
% 

elseif whichone == 3 
ifhildef~=0 

stl = inputCHow many miles will you be driving?'); 
ell = input('At what elevation will you be driving (feet)?'); 
else 

stl = 20; 
ell = 100; 
end 

distvec = 0:.l:stl; 
ifell~=0 

hillvec = [ell*ones(size(l:length(0:.l:stl)))]; 
else 

hillvec = [ 2 eros(size(l:length( 0 :.l:stl)))]; 
end 

type = 'Flat Terrain'; 
hildef=l; 

% or if a terrain profile is to be loaded from a text file. 

% 
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elseif whichone “ 2 

dispCNOTE; The. input file must be set up with two column vectors. The first'); 
dispCcolumn contains the distance from start of simulation in tenths of miles,'); 
dispCand the second column contains the elevation above sea level in feet. The'); 
dispCfile must be a text fiile of the form "filename.dat".'); 
dispC) 

filename = mput('Input the filename: ','s'); 
fid=fopen(filename); 
dat = fscanf(fid,'%g %g',[2inf]); 
fclose(fid); 
distvec= dat(l,:); 
hillvec= dat(2,:); 
type = 'Hilly Terrain'; 
end 

% compute the hill angle (alpha) and output to system 

% 

initelev=hillvec( 1,1); 
for i=2:length(hillvec) 

hillang(i- l)=asin((hillvec(i)-hillvec(i-1))/528); 
end 

hiUang(length(hitlvec))=hillang(length(hillang)); 

% Plot the updated terrain profile for use during simulation 
% 

Figures = get(0,'Children'); 
for i = l:length(Figures) 

if strcmp(get(Figures(i),'Type'), 'figure') 
if strcmp(get(Figures(i), 'Name'), Terrain Profile') 
close(Figures(i)) 
break; 
end 
end 
end 

figureCNumberTitle'.'ofFj'Name'jTerrain Profile',... 

'NextPlot','add’,'Position',[0 0 450 300]); 
plot(distvec,hitlvec),grid,ylabel('feet above sea level') 

axis([0 max(distvec) (min(hillvec)-(.l*min(hillvec))) (max(hillvec)+(.l*max(hitlvec)))]) 
xlabel('distance in miles'),title(type) 

D. FUNCTION EFFFUNC.M TO EXPAND SIZE OF EFFICIENCY MATRIX 
% This function will expand the efficiency matrix in order 
% to decrease interpolation errors. 

% 
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function [efficiency,rpmin,Tlin,nmot,isnew] = efffunc(N,nmot); 

% Efficiency for the 11.3 Hp motor 
% 


isnew= 127; 
if nmot == 1 

isnew=input('Is this the 127 or 86 Unique motor? ’); 
end 


if isnew = 86 

effmat=[35 55 67 68 69 70 71 70 63 55; 

58 68 72 74 75 74 72 70 68 66; 

68 74 78 78.5 78.6 76 74 72 70 69; 

72 77 82 82 81.9 81 80 77 73 70; 

74 79 84 84.4 85 84.8 83 82 80 78; 

78 80 86 87.3 87.2 87.1 86.5 84.5 83 81.5; 

81 84 89 89.7 90.1 89.8 88 87 86 85.5; 

82 86 90.5 92 92.2 92.8 91 89.2 88.5 87; 

83 87 91.7 94 93.5 92.5 91.6 89.8 88.6 87.1; 

84 87.5 93 94.5 94.2 93 91.9 90 89.5 88.5]; 

elseif isnew == 127 

effmat=[10 20 32 39 38 40 51 41 28 18; 

18 36 60 71 70 69 72 63 55 45; 

20 50 76 85.5 82 80.5 79 75 71 64; 

30 60 85.5 86 85.5 84.8 83 82 79 74; 

33 71 86 87.5 86.5 86.5 85.5 85 83.5 81; 

38 80 87 88.5 89 88 87 86 95.5 84; 

40 84 89 90 90.5 90 89 88 87 86; 

50 85 89.5 91 92 91.5 91 90 90.5 89.5; 
50 86 90 92 93 93.5 91 91.5 90.5 89.5; 
48 86.5 90.5 92.5 94 93.5 93 92 91 90.5]; 
nmot= 1; 
end 


[n m]=size(ef&nat); 


% Increases the size of the matrix from n x m to N*n x N*m 
% with linspace interpolation. 

% 

forj=l:n 
for i=2:m 

eftemp(i-l,:)=linspace(effmat(j,i-l),effmat(j,i),N); 
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end 

[c d] = size(eftemp); 
for a=l :c 
for b=l:d-l 

neweffmat(j,(((N-l)*(a-l))+b)) = eftemp(a,b); 
end 
end 
end 


new2=neweffmat'; 
[n m]=size(new2); 


forj=l:n 
for i=2:m 

eftemp2(i-1 ,:)=linspace(new2(j,i- l),new2(j,i),N)5 
end 

[c d] = size(eftemp2); 
for a=l :c 
forb=l:d-l 

neweffmat2(j,(((N-l)*(a-l))+b)) = eftemp2(a,b); 
end 
end 
end 

^ sH***********************************^^******^^**************** 

efficiency = neweffmat2'; 

[n m]=size(efficiency); 

% Corresponding torque and rpm vectors. 

% 

ifisnew==127 % Unique DR127/CR10 
rpmin=linspace(25,5(X)0,n); 

Tlin=linspace(. l,20.25,m); 
elseif isnew = 86 % Unique DR086/CR10 
rpmin=lmspace(200,5600,n); 

Tlin=linspace(.4,15.2,m); 
end 


E. FUNCTION EFFPLTF.M FOR SURFACE PLOT OF EFFICIENCY MATRIX 

% This function performs a surface plot of the motor 
% efficiency 

function effpltf(rpmin,Tlin,effmat,isnew); 
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figure('NumberTide','ofF,'Name','Efficiency Window') 
colormap(jet) 

surf(Tlin,rpmin,effmat),grid 

axis([0 max(Tlin) 0 max(rpmin) 0 100]) 

shading flat 

colorbar 

xlabel('Torque'), 

ylabel('RPM'), 

zlabel('Eff(%)'), 

if isnew =127 

titleCMotor Efficiency for Unique DR127/CR10'); 
elseif isnew == 86 

titleCMotor Efficiency for Unique DR086/CR1 O’); 
end 


F. FUNCTION SPEVAN2.M FOR ON-LINE DISPLAY WINDOW 

function [sys,x0]=spevan2(t,x,u,flag,fignamel); 

% spevan2 S-fUnction for animating the on-line display of dynamic 
% output. To see the specific input corresponding to the parameter 
% u[], open the "On-Line Display" block in the simulation window. 

% Steven Roerig. 7/25/9 

global HNDLl 
if abs(flag) = 2 

if (x(l) ~= Inf) % screen has been initialized 
hndls=get(x(13),'UserData’); 

if u(5) > -.01 & u(5) < .01 % set output to zero 

set(hndls( 1),'String',num2str(0.00)); 
else 

set(hndls( 1) ,'String',num2str(u(5))); 
end 


set(hndls(2),’String',num2str(u(3))); 

set(hndls(3),'String',num2str(u(4))); 

if u(l) > -.001 & u(l) < .001 % set output to zero 

set(hndls(4),'String',num2str(0.00)); 
else 

set(hndls(4),'String',num2str(u( 1))); 
end 

set(hndls(5),'String',num2str(u(2))); 
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set(hndls(8),'Strmg',num2str(u(7))); 
set(hndls(12),’String',num2str(u(6))); 
set(hndls(9),'String',num2str(u(8))); 

if u(9) > -.001 & u(9) < .001 % set vehicle "stopped" 
set(hndls( 10),'S tring',num2str(0.00)); 
set(hndls(ll),'String','Stopped.','Color',[1 1 0]); 
elseif u(9) < -.001 % set vehcile "backward" 

set(hndls(10),'String',num2str(u(9))); 
set(hndls(ll),'String','Going Backward! !','Color',[l 0 0]) 
else % set vehicle "going forward" 

set(hndls(10),'String',num2str(u(9))); 
set(hndls(ll),'String','GoingForward.','Color',[0 1 1]) 

end 

if u(5) > .001 % set drawing current 

set(hndls(6),'XData',l,'YData',2.5,'Color',[0 1 0]); 
set(hndls(7),'String','using'); 
elseif u(5) < -.001 % set regenerating current 

set(hndls(6),'XData',l,'YData',2.5,'Color',[l 0 0]); 
set(hndls(7),'String','regen'); 
else % set motor stopped 

set(hndls(6),'XData',l,'YData',2.5,'Color',[0 0 0]); 
set(hndls(7),'String',''); 
end 

set(hndls(13),'String',num2str(u(10))); 

if u(l 1) <= .01 % set output to zero 

set(hndls(14),'String',num2str(0.00)); 
else 

set(hndls(14),'String',num2str(u(ll))); 

end 

% output the speed to analog display 
y2=[1.5 1.5+1. l*sin((2*pi*u(9)/60))]; 
x2=[10.5 10.5-l.l*cos((2*pi*u(9)/60))]; 
set(hndls(15),'XData',x2,'YData’,y2); 
set(hndls(16),'XData',10.5,'YData',1.5); 
set(hndls(17),'String',num2str(u(12))); 
drawnow; 
end 

sys=[u;x(13)]; % initialize system parameters 

elseif flag = 0 % Initialize the figure "display" window 
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sys=[0;12+l;0;12;0;0]; 

Figures = get(0,'Chil'); 
yesfig = 0; 

for i = l:length(Figures) 
if strcinp(get(Figures(i), 'Name'), 'Output') 
yesfig = 1; 
whichone = i; 
end 
end 

if yesfig == 1 

curhndl = Figures(whichone); 
x0=[lnf;0;0;0;0;0;0;0;0;0;0;0;curhndl]; 
set(curhndl,'Userdata',HNDLl); 
else 

curhndl = animl(fignamel); % draw the speed dial 
cir = pi:-.035:-pi; 
for i=l:length(cir) 
xl(i) = 10.5 + 1.2*cos(cir(i)); 
yl(i) = 1,5 + 1.2*sin(cir(i)); 
end 

for i= 1 :length(cir)/15 
xla(i) = xl(i*length(cir)/12); 
yla(i) = yl(i*length(cir)/12); 
end 

x2=[10.5 9.4]; 
y2=[1.5 1.5]; 

hold on % everything below initializes the display window 
plot(x 1 ,y 1, 'r') ,plot(x 1 a,y 1 a,', w', 'MarkerSize', 11); 
axis([0120 3.8]); 
axis off 

text(8.9,1.55,'0','FontName’,Times New Roman','FontSize',[10]) 
text(10.3,3.,'15','FontName','Times New Roman’,'FontSize',[10]) 
text(l 1.9,1.55,'30','FontName','Times New Roman','FontSize',[10]) 
text(10.4,0,'45','FontName','Times New Roman’,'FontSize’,[10]) 
text(10.1,3.5,'Mph’,'FontName',Times New Roman',’FontSize',[10]) 
text(3.35,3,’Array Current (amps):','FontName',Times New Roman','FontSize',[10]) 
text(3.95,2.5,'Run time (secs):','FontName',Times New Roman','FontSize’,[10]) 
text(3.72,2,'Odometer (miles):','FontName',Times New Roman','FontSize',[10]) 
text(3.5,1.5,'Speedometer (mph);','FontName','Times New Roman','FontSize',[10]); 
text(-.l,3,' Motor Current (amps):','FontName',Times New Roman','FontSize',[10]) 
text(.92,2,'Ahrs Used:','FontName’,'Times NewRoman','FontSize',[10]) 
text(l,1.5,'SOC (%):','FontName’,'Times New Roman’,'FontSize',[10]) 
text(.05,.5,' Vehicle Load (N-m):’,'FontName','Times New Roman','FontSize',[10]); 
text(7.3,3,' RPM;','FontName','Times New Roman','FontSize',[10]); 
text(6.4,2.5,'Hill Grade (%):','FontName','Times New Roman','FontSize',[10]); 



text(-. 1,1,'Battery Voltage (volts):','FontName','Times New Roman’,'FontSize',[10]); 
text(3.35,.5,'Motor Efficiency (%):','FontName','Times New Roman','FontSize',[10]); 
HNDLl(l)=text(2.2,3,' ','EraseMode','background','FontName','Times New 
Roman’,'FontSize',[ 10]); 

HNDLl(2)=text(2.2,2,' ','EraseMode','background','FontName',’Times New 
Roman','FontSize',[ 10]); 

HNDL1 (3)=text(2.2,1.5,' ','EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDL 1 (4)=text(2.2,.5,' ','EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(5)=text(5.6,3,' ','EraseMode’,'background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(6)=plot(l,2.5,'.','MarkerSize',30,'EraseMode','background'); 

HNDL 1 (7)=text( 1.2,2.5,' ’,'EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(8)=text(5.6,2.5,' ',’EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(9)=text(5.6,2,' ','EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(10)=text(5.6,1.5,' ','EraseMode','background','FontName','Times New 
Roman’,'FontSize’,[ 10]); 

HNDLl(ll)=text(4,l,' ','Color',[l 0 0],'EraseMode','background','FontName','Times 
New Roman','FontSize',[10]); 

HNDLl(12)=text(2.2,l,' ’,'EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(13)=text(5.6,.5,' ','EraseMode','background',’FontName',Times New 
Roman',’FontSize',[ 10]); 

HNDLl(14)=text(8,3,' ','EraseMode','background','FontName','Times New 
Roman','FontSize',[ 10]); 

HNDLl(15)=plot(x2,y2,'c','LineWidth',1.5,'EraseMode','background'); 
HNDLl(16)=plot(10.5,1.5,'.','MarkerSize',20,'EraseMode','background'); 
HNDLl(17)=text(8,2.5,' ’,'EraseMode','background','FontName',Times New 
Roman','FontSize',[ 10]); 
hold off 

x0=[lnf;0;0;0;0;0;0;0;0;0;0;0;curhndl]; 

set(curhndl,'UserData',HNDLl); 

end 

else 

sys=n; 

end; 


G. FUNCTION ANIM1.M TO INITIALIZE ON-LINE DISPLAY WINDOW 

function [curhndl] = animl(fignamel); 

%animl Initializes a figure for current animations. 
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% Steven Roerig. 25 July 1994, for SPEV 


% Now initialize the whole figure... 
position=[5 5 825 250]; 
curhndl=figure(... 
'Name',fignamel,... 
'NumberTitle’.'ofF,... 
’BackingStore','off,... 
Position’,position); 
axes(... 

'Units','normalized',... 
'Position',[.05 .1.9 .75],... 
'Visible’,'on',... 
DrawMode','fast'); 


cla reset 
axis off 


H. FUNCTION PLOTS.M TO GENERATE FINAL OUTPUT PLOTS 

% This routine generates the plots for specfice outputs as 
% well as averages for the entire simulation. 

% Plot the miles per hour and the total force on the vehicle 
% 

figure(Position’,[0 340 300 300]) 
subplot(2,1,1 ),plot(miles,mph),grid 
title('miles per hour') 

subplot(2, l,2),plot(miles,Ftot),grid,xlabel('miles driven') 
title('Total force on the vehicle (N)') 

% Plot the current drawn by the motor and the total amp hours used 
% 

figure('Position',[0 0 300 300]) 
subplot(2,l,l),plot(mUes,iq),grid, 
title('current drawn by motor') 

subplot(2,1,2),plot(miles,Ahours),grid,xlabel('miles driven') 
title('Amp-hours used') 

% Plot the power delivered to the drive wheel from the motor and 
% the power required at the rear wheel as a function of the forces 
% 

figure('Position',[340 0 300 300]) 
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subplot(2,1, l),plot(miles,pwrdel),grid, 
title('Power delivered to the rear wheels') 

subplot(2,l,2),plot(miles,pwrreq./((datamat(l,15)/100)*effout)),grid,xlabel('miles driven') 
title('Power required at the rear wheels') 

% Plot the efficiency of the motor for the simulation 
% 

figum('Position',[340 340 400 300]) 
subplot(2,1,1 ),plot(miles ,effout),grid, 
title('Motor efficiency') 

% Display the averages beyond the transient system startup 
% 

subplot(2,l,2) 
axis([0 2 0 2]),axis off 
text(.22,1.5,'Average Speed: 

','FontSize',[9]),text(.8,1.5,num2str(mean(mph(100:length(mph)))),'FontSize',[9]) 
text( 1,1.5 ,'mph','FontSize',[9]) 
text(0,l,'Average Current Drawn: 

','FontSize',[9]),text(.8,l,num2str(mean(iq(100:length(iq)))),'FontSize',[9]) 
text( 1,1 ,'Amps','FontSize',[9]) 
text(-.05,.5,'Average Motor Efficiency: 

’,'FontSize',[9]),text(.8,.5,num2str(mean(effout(100:length(effout)))),'FontSize’,[9]) 

text(l,.5,'%','FontSize',[9]) 

text(1.3,1.5,'Average Load: 

’,'FontSize',[9]),text(l .8,1.5aium2str(mean(n2(100:length(n2))));FontSize',[9]) 

text(2,1.5,'Nm’,'FontSize',[9]) 

text(1.34,l,'Miles Driven: 

’,'FontSize',[9]),text(1.8,l,num2str(miles(length(miles))),'FontSize',[9]) 
text(2,1 ,'miles','FontSize', [9]) 
text(1.25,.5,'Total Ahrs used: 

','FontSize',[9]),text(1.8,.5,num2str(mean(Ahours(100:length(Ahours)))),'FontSize',[9]) 

text(2,.5,'Ahrs',’FontSize',[9]) 
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